Sarcoma microenvironment cell states and ecosystems are associated with prognosis and predict response to immunotherapy.

Characterization of the diverse malignant and stromal cell states that make up soft tissue sarcomas and their correlation with patient outcomes has proven difficult using fixed clinical specimens. Here, we employed EcoTyper, a machine-learning framework, to identify the fundamental cell states and cellular ecosystems that make up sarcomas on a large scale using bulk transcriptomes with clinical annotations. We identified and validated 23 sarcoma-specific, transcriptionally defined cell states, many of which were highly prognostic of patient outcomes across independent datasets. We discovered three conserved cellular communities or ecotypes associated with underlying genomic alterations and distinct clinical outcomes. We show that one ecotype defined by tumor-associated macrophages and epithelial-like malignant cells predicts response to immune-checkpoint inhibition but not chemotherapy and validate our findings in an independent cohort. Our results may enable identification of patients with soft tissue sarcomas who could benefit from immunotherapy and help develop new therapeutic strategies.

Therapy-Induced Senescence Contributes to the Efficacy of Abemaciclib in Patients with Dedifferentiated Liposarcoma.

PURPOSE: We conducted research on CDK4/6 inhibitors (CDK4/6i) simultaneously in the preclinical and clinical spaces to gain a deeper understanding of how senescence influences tumor growth in humans. PATIENTS AND METHODS: We coordinated a first-in-kind phase II clinical trial of the CDK4/6i abemaciclib for patients with progressive dedifferentiated liposarcoma (DDLS) with cellular studies interrogating the molecular basis of geroconversion. RESULTS: Thirty patients with progressing DDLS enrolled and were treated with 200 mg of abemaciclib twice daily. The median progression-free survival was 33 weeks at the time of the data lock, with 23 of 30 progression-free at 12 weeks (76.7%, two-sided 95% CI, 57.7%-90.1%). No new safety signals were identified. Concurrent preclinical work in liposarcoma cell lines identified ANGPTL4 as a necessary late regulator of geroconversion, the pathway from reversible cell-cycle exit to a stably arrested inflammation-provoking senescent cell. Using this insight, we were able to identify patients in which abemaciclib induced tumor cell senescence. Senescence correlated with increased leukocyte infiltration, primarily CD4-positive cells, within a month of therapy. However, those individuals with both senescence and increased TILs were also more likely to acquire resistance later in therapy. These suggest that combining senolytics with abemaciclib in a subset of patients may improve the duration of response. CONCLUSIONS: Abemaciclib was well tolerated and showed promising activity in DDLS. The discovery of ANGPTL4 as a late regulator of geroconversion helped to define how CDK4/6i-induced cellular senescence modulates the immune tumor microenvironment and contributes to both positive and negative clinical outcomes. See related commentary by Weiss et al., p. 649.

Epigenetic age acceleration in surviving versus deceased COVID-19 patients with acute respiratory distress syndrome following hospitalization.

BACKGROUND: Aging has been reported as a major risk factor for severe symptoms and higher mortality rates in COVID-19 patients. Molecular hallmarks such as epigenetic alterations and telomere attenuation reflect the biological process of aging. Epigenetic clocks have been shown to be valuable tools for measuring biological age in various tissues and samples. As such, these epigenetic clocks can determine accelerated biological aging and time-to-mortality across various tissues. Previous reports have shown accelerated biological aging and telomere attrition acceleration following SARS-CoV-2 infection. However, the effect of accelerated epigenetic aging on outcome (death/recovery) in COVID-19 patients with acute respiratory distress syndrome (ARDS) has not been well investigated. RESULTS: In this study, we measured DNA methylation age and telomere attrition in 87 severe COVID-19 cases with ARDS under mechanical ventilation. Furthermore, we compared dynamic changes in epigenetic aging across multiple time points until recovery or death. Epigenetic age was measured using the Horvath, Hannum, DNAm skin and blood, GrimAge, and PhenoAge clocks, whereas telomere length was calculated using the surrogate marker DNAmTL. Our analysis revealed significant accelerated epigenetic aging but no telomere attrition acceleration in severe COVID-19 cases. In addition, we observed epigenetic age deceleration at inclusion versus end of follow-up in recovered but not in deceased COVID-19 cases using certain clocks. When comparing dynamic changes in epigenetic age acceleration (EAA), we detected higher EAA using both the Horvath and PhenoAge clocks in deceased versus recovered patients. The DNAmTL measurements revealed telomere attrition acceleration in deceased COVID-19 patients between inclusion and end of follow-up and a significant change in dynamic telomere attrition acceleration when comparing patients who recovered versus those who died. CONCLUSIONS: EAA and telomere attrition acceleration were associated with treatment outcomes in hospitalized COVID-19 patients with ARDS. A better understanding of the long-term effects of EAA in COVID-19 patients and how they might contribute to long COVID symptoms in recovered individuals is urgently needed.

Immune-related Adverse Events after Immune Checkpoint Blockade-based Therapy Are Associated with Improved Survival in Advanced Sarcomas.

The association between immune-related AEs (irAE) and outcome in patients with sarcoma is not known. We retrospectively reviewed a cohort of patients with advanced sarcoma treated with immune checkpoint blockade (ICB)-based therapy. Association of irAEs with survival was assessed using a Cox regression model that incorporated irAE occurrence as a time-dependent covariate. Tumor samples with available RNA sequencing data were stratified by presence of an irAE to identify patterns of differential gene expression. A total of 131 patients were included. Forty-two (32%) had at least one irAE of any grade and 16 (12%) had at least one grade ≥ 3 irAE. The most common irAEs were hypothyroidism (8.3%), arthralgias (5.3%), pneumonitis (4.6%), allergic reaction (3.8%), and elevated transaminases (3.8%). Median progression-free survival (PFS) and overall survival (OS) from the time of study entry were 11.4 [95% confidence interval (CI), 10.7-15.0) and 74.6 weeks (CI, 44.9-89.7), respectively. On Cox analysis adjusting for clinical covariates that were significant in the univariate setting, the HR for an irAE (HR, 0.662; CI, 0.421-1.041) approached, but did not reach statistical significance for PFS (P = 0.074). Patients had a significantly lower HR for OS (HR, 0.443; CI, 0.246-0.798; P = 0.007) compared with those without or before an irAE. Gene expression profiling on baseline tumor samples found that patients who had an irAE had higher numbers of tumor-infiltrating dendritic cells, CD8+ T cells, and regulatory T cells as well as upregulation of immune and inflammatory pathways. SIGNIFICANCE: irAE after ICB therapy was associated with an improved OS; it also approached statistical significance for improved PFS. Patients who had an irAE were more likely to have an inflamed tumor microenvironment at baseline.

DNA methylation predicts the outcome of COVID-19 patients with acute respiratory distress syndrome.

BACKGROUND: COVID-19 infections could be complicated by acute respiratory distress syndrome (ARDS), increasing mortality risk. We sought to assess the methylome of peripheral blood mononuclear cells in COVID-19 with ARDS. METHODS: We recruited 100 COVID-19 patients with ARDS under mechanical ventilation and 33 non-COVID-19 controls between April and July 2020. COVID-19 patients were followed at four time points for 60 days. DNA methylation and immune cell populations were measured at each time point. A multivariate cox proportional risk regression analysis was conducted to identify predictive signatures according to survival. RESULTS: The comparison of COVID-19 to controls at inclusion revealed the presence of a 14.4% difference in promoter-associated CpGs in genes that control immune-related pathways such as interferon-gamma and interferon-alpha responses. On day 60, 24% of patients died. The inter-comparison of baseline DNA methylation to the last recorded time point in both COVID-19 groups or the intra-comparison between inclusion and the end of follow-up in every group showed that most changes occurred as the disease progressed, mainly in the AIM gene, which is associated with an intensified immune response in those who recovered. The multivariate Cox proportional risk regression analysis showed that higher methylation of the "Apoptotic execution Pathway" genes (ROC1, ZNF789, and H1F0) at inclusion increases mortality risk by over twofold. CONCLUSION: We observed an epigenetic signature of immune-related genes in COVID-19 patients with ARDS. Further, Hypermethylation of the apoptotic execution pathway genes predicts the outcome. TRIAL REGISTRATION: IMRPOVIE study, NCT04473131.

Clinical, genomic, and transcriptomic correlates of response to immune checkpoint blockade-based therapy in a cohort of patients with angiosarcoma treated at a single center.

BACKGROUND: Angiosarcoma is a histologically and molecularly heterogeneous vascular neoplasm with aggressive clinical behavior. Emerging data suggests that immune checkpoint blockade (ICB) is efficacious against some angiosarcomas, particularly cutaneous angiosarcoma of the head and neck (CHN). METHODS: Patients with histologically confirmed angiosarcoma treated with ICB-based therapy at a comprehensive cancer center were retrospectively identified. Clinical characteristics and the results of targeted exome sequencing, transcriptome sequencing, and immunohistochemistry analyses were examined for correlation with clinical benefit. Durable clinical benefit was defined as a progression-free survival (PFS) of ≥16 weeks. RESULTS: For the 35 patients included in the analyses, median PFS and median overall survival (OS) from the time of first ICB-based treatment were 11.9 (95% CI 7.4 to 31.9) and 42.5 (95% CI 19.6 to 114.2) weeks, respectively. Thirteen patients (37%) had PFS ≥16 weeks. Clinical factors associated with longer PFS and longer OS in multivariate analyses were ICB plus other therapy regimens, CHN disease, and white race. Three of 10 patients with CHN angiosarcoma evaluable for tumor mutational burden (TMB) had a TMB ≥10. Five of six patients with CHN angiosarcoma evaluable for mutational signature analysis had a dominant mutational signature associated with ultraviolet (UV) light. No individual gene or genomic pathway was significantly associated with PFS or OS; neither were TMB or UV signature status. Analyses of whole transcriptomes from nine patient tumor samples found upregulation of angiogenesis, inflammatory response, and KRAS signaling pathways, among others, in patients with PFS ≥16 weeks, as well as higher levels of cytotoxic T cells, dendritic cells, and natural killer cells. Patients with PFS <16 weeks had higher numbers of cancer-associated fibroblasts. Immunohistochemistry findings for 12 patients with baseline samples available suggest that neither PD-L1 expression nor presence of tumor-infiltrating lymphocytes at baseline appears necessary for a response to ICB-based therapy. CONCLUSIONS: ICB-based therapy benefits only a subset of angiosarcoma patients. Patients with CHN angiosarcoma are more likely to have PFS ≥16 weeks, a dominant UV mutational signature, and higher TMB than angiosarcomas arising from other primary sites. However, clinical benefit was seen in other angiosarcomas also and was not restricted to tumors with a high TMB, a dominant UV signature, PD-L1 expression, or presence of tumor infiltrating lymphocytes at baseline.

Long-term Follow-up and Patterns of Response, Progression, and Hyperprogression in Patients after PD-1 Blockade in Advanced Sarcoma.

PURPOSE: Programmed cell death protein 1 (PD-1) blockade can mediate objective responses in advanced sarcomas, but their durability has not been established and it is unclear if hyperprogressive disease (HPD) occurs in sarcomas treated with PD-1 inhibitors. EXPERIMENTAL DESIGN: We pooled patients who were treated prospectively with nivolumab or pembrolizumab as monotherapy or with bempegaldesleukin, epacadostat, ipilimumab, or talimogene laherparepvec. We did a new independent assessment for HPD and analyzed clinical, pathologic, and genomic data from baseline tumor biopsies. Our primary endpoint was the incidence of HPD; secondary endpoints were clinical or genomic correlates of response or HPD. RESULTS: We treated 134 patients with advanced sarcoma from 2015 to 2019. Twenty-one patients (16%) had a complete or partial response (CR/PR), and 30% of responses were durable for over 2 years. Forty-eight (36%) patients had stable disease (SD), 45 (34%) had progressive disease without HPD (PD), and 15 (11%) had HPD. Five patients (4%) were not evaluable for HPD. The sarcoma subtypes, sites of metastasis, clinical course, and genomic alterations in patients with PD and HPD were similar, except HPD tumors were smaller at baseline. CONCLUSIONS: In patients with advanced sarcoma, PD-1 blockade can mediate durable responses. HPD occurs in sarcoma at an incidence that is similar to what has been reported in other solid tumors, but patients with HPD were clinically and biologically similar to those who had PD. Further research is required to establish whether HPD is a biologically distinct phenomenon and whether a theoretical risk of HPD should influence patient management.

The derepression of transposable elements in lung cells is associated with the inflammatory response and gene activation in idiopathic pulmonary fibrosis.

BACKGROUND: Transposable elements (TEs) are repetitive sequences of viral origin that compose almost half of the human genome. These elements are tightly controlled within cells, and if activated, they can cause changes in both gene regulation and immune viral responses that have been associated with several chronic inflammatory diseases in humans. As oxidants are potent activators of TEs, and because oxidative injury is a major risk factor in relation to idiopathic pulmonary fibrosis (IPF), we hypothesized that TEs might be involved in the regulation of gene expression and so contribute to inflammation in cases of IPF. IPF is a fatal lung disease that involves the gradual replacement of the alveolar tissue with fibrotic scars as well as the accumulation of inflammatory cells in the lower respiratory tract. Although IPF is known to occur as a result of the complex interaction between age, environmental risk factors (i.e., oxidative stress) and genetics, the relative contributions of these factors to the disease remain unclear. To determine whether TEs are associated with IPF, we compared the transcriptional profiles of the genes and TEs of lung cells obtained from both healthy donors and IPF patients. RESULTS: We quantified TE and gene expression levels using a published bulk RNA-seq dataset containing 24 subjects (16 donors and eight IPF patients), including three lung-cell types per subject, as well as an scRNA-seq dataset concerning 16 subjects (eight donors and eight IPF patients). We found evidence of TE dysregulation in the alveolar type II lung cells and alveolar macrophages of the IPF patients. In addition, the activation of the LINE1 family of elements in IPF is associated with the increased expression of TE cellular regulators (MOV10, IFI16, SAMHD1, and APOBECG3), interferon-stimulating genes (ISG15, IFI6, IFI27, IFI44, and OAS1), chemokines (CX3CL1 and CXCL9), and interleukins (IL15RA). We also propose that TE derepression might be involved in the regulation of previously reported IPF candidate genes (MUC5B, CHL1, SPP1, and MMP7). CONCLUSION: Based on our findings, we propose that TE derepression plays an important role in the regulation of gene expression and can also prompt both the recruitment of inflammatory processes and the disruption of the immunological balance, which can lead to chronic inflammation in IPF.

Whole-methylome analysis of circulating monocytes in acute diabetic Charcot foot reveals differentially methylated genes involved in the formation of osteoclasts.

AIM: To assess whether DNA methylation of monocytes play a role in the development of acute diabetic Charcot foot (CF). PATIENTS & METHODS: We studied the whole methylome (WM) of circulating monocytes in 18 patients with Type 2 diabetes (T2D) and acute CF, 18 T2D patients with equivalent neuropathy and 18 T2D patients without neuropathy, using the enhanced reduced representation bisulfite sequencing technique. RESULTS & CONCLUSION: WM analysis demonstrated that CF monocytes are differentially methylated compared with non-CF monocytes, in both CpG-site and gene-mapped analysis approaches. Among the methylated genes, several are involved in the migration process during monocyte differentiation into osteoclasts or are indirectly involved through the regulation of inflammatory pathways. Finally, we demonstrated an association between methylation and gene expression in cis- and trans-association.

Single-dose versus 7-day-dose metronidazole for the treatment of trichomoniasis in women: an open-label, randomised controlled trial.

BACKGROUND: Among women, trichomoniasis is the most common non-viral sexually transmitted infection worldwide, and is associated with serious reproductive morbidity, poor birth outcomes, and amplified HIV transmission. Single-dose metronidazole is the first-line treatment for trichomoniasis. However, bacterial vaginosis can alter treatment efficacy in HIV-infected women, and single-dose metronidazole treatment might not always clear infection. We compared single-dose metronidazole with a 7-day dose for the treatment of trichomoniasis among HIV-uninfected, non-pregnant women and tested whether efficacy was modified by bacterial vaginosis. METHODS: In this multicentre, open-label, randomised controlled trial, participants were recruited at three sexual health clinics in the USA. We included women positive for Trichomonas vaginalis infection according to clinical screening. Participants were randomly assigned (1:1) to receive either a single dose of 2 g of metronidazole (single-dose group) or 500 mg of metronidazole twice daily for 7 days (7-day-dose group). The randomisation was done by blocks of four or six for each site. Patients and investigators were aware of treatment assignment. The primary outcome was T vaginalis infection by intention to treat, at test-of-cure 4 weeks after completion of treatment. The analysis of the primary outcome per nucleic acid amplification test or culture was also stratified by bacterial vaginosis status. This trial is registered with ClinicalTrials.gov, number NCT01018095, and with the US Food and Drug Administration, number IND118276, and is closed to accrual. FINDINGS: Participants were recruited from Oct 6, 2014, to April 26, 2017. Of the 1028 patients assessed for eligibility, 623 women were randomly assigned to treatment groups (311 women in the single-dose group and 312 women in the 7-day-dose group; intention-to-treat population). Although planned enrolment had been 1664 women, the study was stopped early because of funding limitations. Patients in the 7-day-dose group were less likely to be T vaginalis positive at test-of-cure than those in the single-dose group (34 [11%] of 312 vs 58 [19%] of 311, relative risk 0·55, 95% CI 0·34-0·70; p<0·0001). Bacterial vaginosis status had no significant effect on relative risk (p=0·17). Self-reported adherence was 96% in the 7-day-dose group and 99% in the single-dose group. Side-effects were similar by group; the most common side-effect was nausea (124 [23%]), followed by headache (38 [7%]) and vomiting (19 [4%]). INTERPRETATION: The 7-day-dose metronidazole should be the preferred treatment for trichomoniasis among women. FUNDING: National Institutes of Health.

Genetic Indicators of Drug Resistance in the Highly Repetitive Genome of Trichomonas vaginalis.

Trichomonas vaginalis, the most common nonviral sexually transmitted parasite, causes ∼283 million trichomoniasis infections annually and is associated with pregnancy complications and increased risk of HIV-1 acquisition. The antimicrobial drug metronidazole is used for treatment, but in a fraction of clinical cases, the parasites can become resistant to this drug. We undertook sequencing of multiple clinical isolates and lab derived lines to identify genetic markers and mechanisms of metronidazole resistance. Reduced representation genome sequencing of ∼100 T. vaginalis clinical isolates identified 3,923 SNP markers and presence of a bipartite population structure. Linkage disequilibrium was found to decay rapidly, suggesting genome-wide recombination and the feasibility of genetic association studies in the parasite. We identified 72 SNPs associated with metronidazole resistance, and a comparison of SNPs within several lab-derived resistant lines revealed an overlap with the clinically resistant isolates. We identified SNPs in genes for which no function has yet been assigned, as well as in functionally-characterized genes relevant to drug resistance (e.g., pyruvate:ferredoxin oxidoreductase). Transcription profiles of resistant strains showed common changes in genes involved in drug activation (e.g., flavin reductase), accumulation (e.g., multidrug resistance pump), and detoxification (e.g., nitroreductase). Finally, we identified convergent genetic changes in lab-derived resistant lines of Tritrichomonas foetus, a distantly related species that causes venereal disease in cattle. Shared genetic changes within and between T. vaginalis and Tr. foetus parasites suggest conservation of the pathways through which adaptation has occurred. These findings extend our knowledge of drug resistance in the parasite, providing a panel of markers that can be used as a diagnostic tool.

The impact of genomics on population genetics of parasitic diseases.

Parasites, defined as eukaryotic microbes and parasitic worms that cause global diseases of human and veterinary importance, span many lineages in the eukaryotic Tree of Life. Historically challenging to study due to their complicated life-cycles and association with impoverished settings, their inherent complexities are now being elucidated by genome sequencing. Over the course of the last decade, projects in large sequencing centers, and increasingly frequently in individual research labs, have sequenced dozens of parasite reference genomes and field isolates from patient populations. This 'tsunami' of genomic data is answering questions about parasite genetic diversity, signatures of evolution orchestrated through anti-parasitic drug and host immune pressure, and the characteristics of populations. This brief review focuses on the state of the art of parasitic protist genomics, how the peculiar genomes of parasites are driving creative methods for their sequencing, and the impact that next-generation sequencing is having on our understanding of parasite population genomics and control of the diseases they cause.

The Tc1/mariner transposable element family shapes genetic variation and gene expression in the protist Trichomonas vaginalis.

BACKGROUND: Trichomonas vaginalis is the most prevalent non-viral sexually transmitted parasite. Although the protist is presumed to reproduce asexually, 60% of its haploid genome contains transposable elements (TEs), known contributors to genome variability. The availability of a draft genome sequence and our collection of >200 global isolates of T. vaginalis facilitate the study and analysis of TE population dynamics and their contribution to genomic variability in this protist. RESULTS: We present here a pilot study of a subset of class II Tc1/mariner TEs that belong to the T. vaginalis Tvmar1 family. We report the genetic structure of 19 Tvmar1 loci, their ability to encode a full-length transposase protein, and their insertion frequencies in 94 global isolates from seven regions of the world. While most of the Tvmar1 elements studied exhibited low insertion frequencies, two of the 19 loci (locus 1 and locus 9) show high insertion frequencies of 1.00 and 0.96, respectively. The genetic structuring of the global populations identified by principal component analysis (PCA) of the Tvmar1 loci is in general agreement with published data based on genotyping, showing that Tvmar1 polymorphisms are a robust indicator of T. vaginalis genetic history. Analysis of expression of 22 genes flanking 13 Tvmar1 loci indicated significantly altered expression of six of the genes next to five Tvmar1 insertions, suggesting that the insertions have functional implications for T. vaginalis gene expression. CONCLUSIONS: Our study is the first in T. vaginalis to describe Tvmar1 population dynamics and its contribution to genetic variability of the parasite. We show that a majority of our studied Tvmar1 insertion loci exist at very low frequencies in the global population, and insertions are variable between geographical isolates. In addition, we observe that low frequency insertion is related to reduced or abolished expression of flanking genes. While low insertion frequencies might be expected, we identified two Tvmar1 insertion loci that are fixed across global populations. This observation indicates that Tvmar1 insertion may have differing impacts and fitness costs in the host genome and may play varying roles in the adaptive evolution of T. vaginalis.

The population genomics of repeated evolution in the blind cavefish Astyanax mexicanus.

Distinct populations of Astyanax mexicanus cavefish offer striking examples of repeatable convergence or parallelism in their independent evolutions from surface to cave phenotypes. However, the extent to which the repeatability of evolution occurred at the genetic level remains poorly understood. To address this, we first characterized the genetic diversity of 518 single-nucleotide polymorphisms (SNPs), obtained through RAD tag sequencing and distributed throughout the genome, in seven cave and three groups of surface populations. The cave populations represented two distinct lineages (old and new). Thirty-one SNPs were significantly differentiated between surface and old cave populations, two SNPs were differentiated between surface and new cave populations, and 44 SNPs were significantly differentiated in both old and new cave populations. In addition, we determined whether these SNPs map to the same locations of previously described quantitative trait loci (QTL) between surface and cave populations. A total of 25 differentiated SNPs co-map with several QTL, such as one containing a fibroblast growth factor gene (Fgf8) involved in eye development and lens size. Further, the identity of many SNPs that co-mapped with QTL was the same in independently derived cave populations. These conclusions were further confirmed by haplotype analyses of SNPs within QTL regions. Our findings indicate that the repeatability of evolution at the genetic level is substantial, suggesting that ancestral standing genetic variation significantly contributed to the population genetic variability used in adaptation to the cave environment.

Getting trichy: tools and approaches to interrogating Trichomonas vaginalis in a post-genome world.

Trichomonas vaginalis is a parasite of the urogenital tract in men and women, with a worldwide presence and significant implications for global public health. T. vaginalis research entered the age of genomics with the publication of the first genome sequence in 2007, but subsequent utilization of other 'omics' technologies and methods has been slow. Here, we review some of the tools and approaches available to interrogate T. vaginalis biology, with an emphasis on recent advances and current limitations, and draw attention to areas where further efforts are needed to examine effectively the complex and intriguing biology of the parasite.

Gene flow and population structure in the Mexican blind cavefish complex (Astyanax mexicanus).

BACKGROUND: Cave animals converge evolutionarily on a suite of troglomorphic traits, the best known of which are eyelessness and depigmentation. We studied 11 cave and 10 surface populations of Astyanax mexicanus in order to better understand the evolutionary origins of the cave forms, the basic genetic structuring of both cave and surface populations, and the degree to which present day migration among them affects their genetic divergence. RESULTS: To assess the genetic structure within populations and the relationships among them we genotyped individuals at 26 microsatellite loci. We found that surface populations are similar to one another, despite their relatively large geographic separation, whereas the cave populations are better differentiated. The cave populations we studied span the full range of the cave forms in three separate geographic regions and have at least five separate evolutionary origins. Cave populations had lower genetic diversity than surface populations, correlated with their smaller effective population sizes, probably the result of food and space limitations. Some of the cave populations receive migrants from the surface and exchange migrants with one another, especially when geographically close. This admixture results in significant heterozygote deficiencies at numerous loci due to Wahlund effects. Cave populations receiving migrants from the surface contain small numbers of individuals that are intermediate in both phenotype and genotype, affirming at least limited gene flow from the surface. CONCLUSIONS: Cave populations of this species are derived from two different surface stocks denoted "old" and "new." The old stock colonized caves at least three times independently while the new stock colonized caves at least twice independently. Thus, the similar cave phenotypes found in these caves are the result of repeated convergences. These phenotypic convergences have occurred in spite of gene flow from surface populations suggesting either strong natural or sexual selection for alleles responsible for the cave phenotype in the cave environment.

Genotyping with Sequenom.

Often in evolutionary genetics research, one needs to analyze polymorphisms in populations for which cost-efficient high-throughput arrays are nonexistent, either because the species is not a model organism or because the populations have been subjected to such specific conditions that their base variation is almost unique. In this situation, custom-made genotyping assays are required. Sequenom's MassARRAY(®) genotyping platform is a powerful and flexible method for assaying up to a few thousand markers and up to thousands of individuals. It is based on distinguishing allele-specific primer extension products by mass spectrometry (MALDI-TOF). Most stages of the experimental protocol reflect adaptations of established PCR protocols to multiplexing, which allows the simultaneous amplification and detection of multiple markers per reaction.

Experimental evolution reveals natural selection on standing genetic variation.

Evolution depends on genetic variation generated by mutation or recombination from standing genetic variation. In sexual organisms, little is known about the molecular population genetics of adaptation and reverse evolution. We carry out 50 generations of experimental reverse evolution in populations of Drosophila melanogaster, previously differentiated by forward evolution, and follow changes in the frequency of SNPs in both arms of the third chromosome. We characterize the effects of sampling finite population sizes and natural selection at the genotype level. We demonstrate that selection has occurred at several loci and further that there is no general loss or gain of allele diversity. We also observe that despite the complete convergence to ancestral levels of adaptation, allele frequencies only show partial return.